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1 Introduction 

Motivated by recent experimental progress to manipulate the refractive index of dielectric ma- 
terials by strong laser beams, we study some aspects of the quantum radiation created by such 
refractive index perturbations. The Kerr effect in non-linear dielectrics describes a change of 
the refractive index n due to electromagnetic radiation with the intensity / 

n = n + n 2 I, (1) 

where uq is the unperturbed index of refraction and n 2 is the Kerr coefficient. For example, 
for fused silica, it is n 2 = 3 x 10 -16 W -1 cm 2 . Since the intensity / will be far below l/n 2 , 
the change 5n in the refractive index n will be small. However, with a space-time dependent 
intensity I(t,r), we can generate a varying refractive index perturbation 

n(t,r) = n + 5n(t,r) , 5n < 1 . (2) 

Such a perturbation is able to create photon pairs out of the vacuum - more precisely, the ground 
state of the quantized electromagnetic field in the dielectric medium. In order to investigate this 
effect, we will make the following simplifications: 

• We assume zero temperature and thus start in the ground state. 

• We omit the photon polarizations and consider scalar fields. 

• We neglect the dispersion n(uj) of the medium. 

• Since 5n <C 1 is small, we use perturbation theory. 

The first two assumptions are not critical and not difficult to fix. However, taking into account 
the dispersion n(u) of the medium is rather non-trivial and calculations beyond perturbation 
theory are also far more involved in general. 



2 Perturbation Theory 



Given the aforementioned assumptions, we may start from the Lagrangian (H = cq = 1) 

C = \ (eE 2 -B 2 ) = \ ([n + 5n] 2 E 2 - B 2 ) , (3) 

where e is the dielectric permittivity of the medium and E and B are the microscopic electric 
and magnetic fields, respectively. Now standard quantization procedure yields the Hamiltonian 

H = \ [ d 3 r( [ D * 2 +B 2 ) =60 + 6,, (4) 
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where D = eE is the electric displacement field including the polarization P of the medium 
D = E + P. Using 5n <C 1, we expand H into powers of 5n such that the lowest order 
Hq is the undisturbed Hamiltonian and the rest H\ contains the effect of the refractive index 
perturbation. Then we may use usual time-dependent perturbation theory in the interaction 
picture to calculate the final quantum state l^out) starting from the initial vacuum state |0) 

|*out> = \0)-iJ dtH x {t) |0) + 0{5n 2 ) . (5) 

Inserting the expansion into creation and annihilation operators for the quantized electromag- 
netic field, we find the production of photon pairs (to lowest order in 5n) 

|*out> = |0> + A k , k > \k, k') + 0(5n 2 ) , (6) 
fc,fc' 

with the two-photon amplitude/probability (see, e.g., pQ and [2]) 



Here k and k' are the wave-numbers of the created photon pair and ujk oj' k their frequencies. 



A k y\=^ 5n{u: k +^k + k') . (7) 
n 



3 Single One-Parameter Pulse 

As a first example, let us consider a simple pulse whose spatial and temporal extent are given 
by the same parameter fi, i.e., 

5n{t,r) = 5nf(nt,nr/c), (8) 

with a well behaved function / = 0(1) (e.g., Gaussian). Then it turns out that the total emission 
probability is independent of 

P = EI" 4 ^'I 2 = °(^ 2 )« 1 ' (9) 
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where the precise pre-factor in front of 5n 2 depends on the precise pulse shape. The typical 
photon energy, however, is given by 

E = ^J2\ A ^'\ 2u k = o(n). (io) 

k,k' 

Nevertheless, the total probability is bound to be small. For an experimental verification, 
therefore, a pulse with more than one parameter might be more useful. 

4 Single Two-Parameter Pulse 

Thus, as a second example, let us assume that the spatial and temporal scales of the pulse are 
given by two different parameters, Cl\ and CI2 

8n(t,r) = 5nf(n 1 t,tt 2 r/c). (11) 

If ili and Q2 are of the same order, we basically reproduce the previous case. Thus, let us 
consider the two limiting cases fii <C ^2 and Cli 3> O2. 

• For CI 1 <C ^2, the spatial extend is very small and we have a "point-like pulse". 

• For $7i 3> the pulse is almost spatially homogeneous and thus the situation is very 
similar to "cosmological particle creation". 

In the first case ("point-like pulse"), one can derive a closed expression for the average energy 
emitted by the refractive index perturbation [1J 

(E) x Jdt(j^y d 3 r Sn(t, r)j ) ' = O (J| 5n^j . (12) 

However, since the typical photon energy is set by the temporal frequency Cli and not by the 
spatial scale fl 2 , the total probability is extremely small 

O 6 



E = OiSlx) ~> P = O ( ^|5n 2 ) < 1. (13) 



Thus this effect is also probably very hard to detect. 

5 Cosmological Particle Creation 

In the opposite limit Q,\ S> 02, we may approximate the pulse with 

Sn(t, r) = 5n /(fitf, fl 2 r/c) (14) 
by a effectively spatially homogeneous profile Sn(t, r) ~ 8n(t) and thus the Lagrangian reads 

C~\(e(t)E 2 -B 2 ) = \(n*(t)E 2 -B 2 ) . (15) 
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This is completely analogous to an expanding/contracting universe where the scale factor a 2 (t) 
corresponds to the refractive index n(t). Hence we can immediately apply the results known for 
this situation, see, e.g., [3] and references therein. 

• Photons are created in pairs with nearly opposite momenta k' ~ —k. 

• Thus, they are in an entangled (squeezed) state. 

• The typical photon energy scales with E = 0{Vt\). 

• The total probability is enlarged by a volume enhancement factor f^/f^ 
In view of this enhancement, this effect might well be observable. 




6 Mixed Case 

For completeness, let us briefly discuss the scenario with three different parameters 

5n(t, r) = 5n f(0,\t, f^/c, ^y/c, Slzz/c) , (17) 

where one spatial scale, say Q2, is much smaller than the temporal frequency fii while the other, 
is much larger - corresponding to a long "needle-like pulse" pointing in x-direction 

fi 2 «^i<^3- (18) 

Again, the typical photon energy is set by the temporal frequency f^i (which is the generic case). 
The total probability is enhanced by a length factor f^/f^ but suppressed by the small cross 
section area oc of the "needle" 

^°(l| Sfi2 )' (19) 

7 Moving Pulse 

So far, we have discussed the situation of a refractive index perturbation 5n(t, r) appearing in 
a localized spatial region and then disappearing again. Now, let us consider the scenario of a 
moving pulse - where we first study the case of constant velocity v 

5n(t,r) = 5nf{Q[r-vt]/c). (20) 

The resulting quantum radiation crucially depends on whether the velocity v is smaller or larger 
than the speed of light c = 1/n in the medium, see also |4]. 
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• For a "sub-luminal" pulse v < c, there is no effect (to lowest order in 5n). 

This can be understood in the following way: electrodynamics in a (homogeneous and 
isotropic) dielectric medium is formally invariant under "Lorentz" transformations with 
the speed of light cq in vacuum being replaced by the speed of light c = 1/n in the 
medium. Thus, after such a "Lorentz" boost, the refractive index perturbation Sn(t, r) 
becomes effectively stationary 5n(t, r) — > 5n{r) and thus - as we have seen before - does 
not produce radiation. 

• For a "super-luminal" pulse v > c, we get "quantum Cherenkov" radiation. 

In this case, a suitable "Lorentz" boost yields an instantaneous perturbation 5n[t) x f(r). 

In the case of a moving pulse (without beginning and end), we cannot compute a total emis- 
sion probability but only the emission probability per unit time (as in Fermi's golden rule, for 
example). In the "super-luminal" case, the situation after the "Lorentz" boost is analogous to 
the previous section with i}\ = O3 — >• and - > and thus we get 

^ = 0{Q5n 2 ). (21) 

Due to phase matching conditions, the photon pairs are predominantly emitted in forward 
direction within an emission angle 7? = 0{\Jv — c). This angle closes for v ^ c and thus the 
emission probability (per unit time) vanishes in this limit - consistent with the above picture. 
As usual, the typical photon energy scales with E = 0(Q). 



8 Hawking Radiation 

Now let us go right to the borderline between the "sub-luminal" and the "super-luminal" case 
discussed above and study a pulse which is slower than the medium speed of light outside the 
pulse c = 1/n but faster than the medium speed of light inside l/(n + Sn) 

5n(t,r) = 6nf(n[r-vt]/c), c = ->v> — 3—. (22) 

n n + on 

In this case, we get the analogue of black hole horizon (light cannot escape) at the front end 
of the pulse and the analogue of a white hole horizon (light cannot enter) at its back and, see 
Fig. [TJ The idea of creating analogues of black (or white) holes in the laboratory and so to 
experimentally test Hawking radiation [6] goes back to Bill Unruh [7], see also [5J [TO] , 
Originally, the analogy was developed for sound waves [7], but later electromagnetic waves were 
considered as well. In [HUE] proposals based on "slow light" (i.e., electromagnetically induced 
transparency) where put forward. However, these suggestions had various problems, see, e.g., 
[131 El ES] . A set-up based on moving dielectrics was proposed in [16], but there the difficulty 
was to have the medium flowing faster than the speed of light in the medium - which is typically 
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Figure 1: Sketch of the pulse creating analogues of black and white hole horizons. 

very fast. Later is was realized [17] that it is not really necessary to actually move the medium - 
instead a moving pulse which changes the local propagation speed can be enough, see also |18j . 
The Hawking temperature (in the co-moving frame) is determined by the gradient of the refrac- 
tive index, see, e.g., [71 [T6] 

^Hawking °c V<5n ~ £1 5n . (23) 

Since Hawking radiation is thermal, the total emission probability (per unit time) scales as 

^ = 0(ATl awking )^0{nSn 3 ) . (24) 

As we may infer from the 5n ?J scaling, this effect is beyond lowest order perturbation theory. 
In one spatial dimension (such as in an optical fibre, see [18]), we would get P/T oc 5n. This 
shows that this effect is non-perturbative - since the perturbative expansion always starts with 
8n 2 . This can be explained by the fact that Hawking radiation occurs due to the tearing apart 
of modes at the black-hole horizon, which takes a time duration long compared to 1/(0 Sn). 
An experiment with a pulse satisfying Eq. (|22p in some frequency range has been done [19] and 
it has been claimed that this was the first observation of the analogue of Hawking radiation. 
However, comparison with the above estimates casts some doubts at these claims, see [201 121j . 



9 Non-Inertial Pulse 

After having discussed the case of constant velocity, let us turn to non-inertial pulse motion 

Sn(t,r) = 5nf(n[r-r P (t)]/c), (25) 

where rp(t) is the trajectory of the pulse. For simplicity, let us consider the case of approximately 
uniform acceleration rp ~ const for some time. In this case, the quantum radiation created by 
the refractive index perturbation can be nicely understood as a signature of the Unruh effect. 
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The Unruh effect [22] describes the striking prediction that a uniformly accelerated detector 
experiences the inertial Minkowski vacuum state as thermal bath with Unruh temperature 

h 

TUnruh = — -: a , (26) 

2irk B c 

where a is the acceleration. Unfortunately, due the the factors h and Co, this temperature is quite 
low for every-day accelerations - e.g., the earth's gravitational acceleration a = g = 9.81 m/s 2 
corresponds to T Unruh = O(10~ 20 K). This is one of the reasons why this effect has not been 
directly observed yet, see also [231 I2H I2S] and references therein. 

thermal bath accelerated scatterer 

v -\ >o ' / / ... Y : o o 

scatterer i u nr uh -^-^ a 

Figure 2: Sketch of the scattering of two photons (red arrows) in the accelerated frame (left) 
and translation of this event to the inertial frame (right). 

However, it should be possible to observe a signature of the Unruh effect with accelerated 
refractive index perturbations, for example. As discussed before, the ground state of a dielectric 
medium behaves analogous to the Minkowski vacuum with the vacuum speed of light cq being 
replaced by the medium light velocity c. Thus, in its non-inertial frame, the accelerated refractive 
index perturbation experiences the ground state of the medium as a thermal bath. Then there 
is a finite probability that one photon out of this thermal bath is scattered by the refractive 
index perturbation into another mode. This scattering event in the accelerated frame, when 
translated back into the inertial (laboratory) frame, corresponds to the emission of real photon 
pairs, as sketched in Fig. [2j In agreement with this intuitive picture, the emission probability 
per unit time scale as (again using perturbation theory) 

P ( 5n 2 \ 

O (cr scat tering?Unruh) ~> ° TTT * P) ' ( 27 ) 



— ^ ^scattering--- unruh,/ \ Q 2 

Note that this expression is valid for 3> |f\p| only since otherwise the internal width ~ of 
the pulse smears out its trajectory rp(t) and thus the picture based on the Unruh effect does 
not apply any-more (even though there still would be quantum radiation). 



10 Dynamical Casimir Effect 

Finally, let us discuss the relation between the quantum radiation given off by these refractive 
index perturbations and the dynamical Casimir effect. In the original setting, the static Casimir 
effect [26] describes the attraction (or repulsion - depending on the boundary conditions) of 
two conducting plates at rest in vacuum, which is caused by the distortion of the quantum 
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vacuum fluctuations of the electromagnetic field. The dynamical Casimir effect then refers to 
the creation of photon pairs out of QED vacuum due to the non-inertial motion of one or both 
mirrors, see, e.g., [27] . 

Replacing the two mirrors by two bodies of dielectric material, one can also get a static Casimir 
attraction. Thus, the quantum radiation created by the non-inertial motion of a dielectric body 
in vacuum could also be called dynamical Casimir effect. As we have discussed before, this 
scenario is formally equivalent to a refractive index perturbation moving in a dielectric medium. 
Thus, the signatures of the Unruh effect caused by the non-inertial motion of such a perturbation 
as discussed in the previous Section can also be viewed as a manifestatiorfj] of the dynamical 
Casimir effect. 
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